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INTRODUCTION 

Although similar visualization materials and techniques have been discussed in previous work by Hendricks 
et al. (1997), we realized that printed frames, grabbed from videotapes of the flows in beds of packed spheres, failed 
to reproduce what was visualized. Herein we describe the events associated with packed beds of twisted tapes, but 
to fully appreciate the complexity of the flow fields, it becomes necessary to watch the videotape recording. How- 
ever, the Full Flow Field Tracking (FFFT) method (Braun et al., 1988) can be applied to visualize and quantize the 
flow patterns and fluid velocities within a packed bed or a porous medium. 

Packed beds of twisted tapes may serve as an alternative to porous-media packed beds in heat pipe applications 
in low- and high-bodv force fields, such as in space and gas turbine applications. Potentially, twisted tapes could 
also function as reaction surfaces where uniformity of mixing is sought, in this study the test section assembly 
simulated a canister of twisted tapes for heat, mass, and reaction exchange. An assembled cylinder would contain 6 
to 10 such canisters in series. In turn, several cylinders would be bundled into an array. We tested only one simu- 
lated canister, although for any practical application simulation of an entire cylinder and array would be necessary. 
Power developed along the path and heat transfer and pressure drops downstream would modify the flows in the 
upstream leading canister even to the point of choking the element. Choking is quite serious because the power 
generation in a practical system is nearly constant and failure of the element becomes imminent. 


ANALYSIS 


Using twisted tapes (fig. 1, from Smithberg and I .and is, 1964) is a well-known method for augmenting heat 
transfer in tubes at the expense of pressure drop in single-phase flows (Hong and Bergles, 1976; Lopina and Bergles, 
1969; Bergles, 1998). Yet at constant pump power, twisted-tape, swirl-flow heat transfer can be increased by 20% 
over that of a straight tube (Lopina and Bergles, 1969). We will use these references as the starting point for 
representing the ideal packed bed of N twisted tapes (cf figs. 1 and 2) and work toward a porous-media model. The 
details are presented in appendix A. 

For a single twisted tape in a tube the tangential fluid velocity is usually assumed to be linear with radial 
position (rotating slug flow) or 


v e = 


27 lUr 
Ho 


( 1 ) 
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v, = ■ 


kUD 

H n 


where r is the radial position, D is the twisted-tape diameter, H 0 is the twist through 360° (one full wave ), and JJ is 
the bulk average axial velocity. (All symbols are defined in appendix B.) 


U = 


w 



( 3 ) 


where W is the mass flow rate, p is the average fluid bulk density, t is the tape thickness, and w is the tape width, and 
for t « D, wt —> Dt. 


Single Twisted Tape in Tube 

From figure 2 of Smithberg and Landis (1964), the data closely follow' equation (1), except at the wall and 
centerline where v 0 = 0, and nearly fit the following form: 


2Uc 
u c 
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( \ 

27" 

H 
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( 4 ) 


which is a close approximation to equation (1). For these data, 1 Re D = 137 000, HJD = 3.62, D = 3.51 cm 
(1.382 in.), and v 0 =71.63 m/s (235 ft/s) at r= 1.65 cm (0.65 in.). At that point r/r 0 = 0.94, v 0 = K(2r/D)(D/H 0 )U, 
or U = 82.3 m/s (270 ft/s). From equation (4) 


C, 


235 

270 

0.763 


1.14 


( 5 ) 


Here v e -> U and even for this case neither an effective velocity (eq. (6) with C 0 ~ C e ) nor an effective flow path 
(eq. (A12b)) is sufficient to account for the measured pressure drop increase in terms of flf 0 . 

For laminar flows with Reynolds numbers less than 150, Date (1974) determined that D!H 0 corrections are not 
required. However, for 150 < Re < 2000 (and for turbulent flows), DIH 0 corrections are necessary. From th eflf 0 data 
of Smithberg and Landis (1964), C e — > 3 with an equivalent velocity defined by using equation (2), 


Ue 

u 



( 6 ) 


The simple empirical form (eq. (6)) tends to group the turbulent friction data of Smithberg and Landis (1964) 
at a higher Reynolds number and identifies ( DIH () ) 2 as a significant parameter for the analysis. However, there are 
additional Reynolds number and surface roughness dependencies (Gambill and Bundy, 1962). For example, at 


'The data point Re D = 137 000 with DID H = 1.716 becomes Re c = 79 837. Extrapolating the HJD = 3.62 locus may give 
f/f g = 2.8, and this ratio may also be affected by roughness. 
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Re = 25 000 the agreement of friction factors is good at high HJD, but at Re = 6000 the agreement is better at low 
HJD but is still not that good. The surface roughness is not given and is assumed to be that of a commercial tube. 

Gambill and Bundy (1962) correlate the isothermal, single-twisted-tape data of several investigators with 
different surface roughness 8 ID g factors. 


where 



(A43) 




8 T 

n = 0.81 exp 

-1700 

, A? J 


(A43a) 


fa=4f 0 = 



(A43b) 


and y is the number of tube diameters per 1 80° of twist (y = HID or 2y = HJD). 


HJD = 2y 

5 ID e = 0 

h!D e = 0.00005 

m e = 0.0005 
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f a = 4 f 0 = 4(0. 046/Re 0 ' 2 ) :f s from equation (A43) 
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Packed Bed of Twisted Tapes 

In the packed-bed experiment described herein the tapes were not bounded by tube walls. They were assembled 
into a uniform matrix of twisted tapes with the same twist direction, twist-to-diameter ratio H!D 0 , and thickness t. 
The boundary conditions changed from those of a single twisted tape. The tangential velocity became zero at the 
center of the twisted tape and at the tangent points of the packed bed of virtual tubes of diameter D 0 because the 
velocity fields were counterrotating. The rotating velocity external to the confines of the virtual cylinders was 
assumed to be small, thus permitting a local region of axial flow. As noted earlier, corrections for H!D 0 effects in 
low-Reynolds-number flows are not required (Date, 1 974). For HID a > 3.6 and w — » D, the tape lengths before and 
after the twist do not differ substantially. 

In packed beds the superficial velocity is related to the bed porosity by 



e 

where U Q is the empty or unpacked bed velocity and e is related to the bed volume V as 

g _ j _ Xsolid 
Hotal 


(7) 


( 8 ) 
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For the present case each of 48 twisted tapes was considered as encased in a virtual tube, where the tape width 
w ----- 0.3234 cm (0. 1273 in.) and the tape thickness t = 0.1275 cm (0.0502 in.). 


D o=\W 



(A4) 


Thus, D 0 = 0.348 cm (0.137 in.) and H Q = 5.503 cm (2.17 in.). For this geometry HJD 0 = 15.8 -> HID Q = 7.9, or 
three 360° twists in 16.5 cm (6.5 in.). 

For the packed bed of 48 twisted tapes considered herein, the orientation of the twist angle was somewhat 
irregular (figs. 2(b), (d), (e), and (f)); however, by using the prior relation for porosity and variations in t and tv, the 
estimated porosity range becomes (see measurements and estimates in appendix A) 

0.59 < e<0.63 (9) 

with a mass flow range of 

1.6 <-^^< 1.7 ( 10 ) 

W 


where A Q is the cross section of the empty tube. 


Parameters for Single Twisted Tape in Tube 


The analysis and data of Smithberg and Landis (1964) were considered to provide insights into the effects of 
twist HID or H/D 0 and pressure drop. Reformulating the Ergun parameter T Ergun (Ergun, 1952) in terms of the 
Fanning friction factor (see appendix A) gives 

•/Smithberg-Landis = ’J^Ergim (A14) 


'Ergun 


'pAP' 

f " 3 1 

1) \ 
u p 

i Go J 

[i-ej 

, L 0 j 


where 


(A16) 


G 0 =pU 0 


At high Reynolds numbers 


Igrgun — ^ Constant — > 0.014 


for e -> 1 and Re » 2000 


(A15) 


(A25) 


where 0.014 is the commercial rough-tube equivalent. In equation (A 16), L g represents the straight-line distance 
between pressure taps. The twisted-tape length L varies little from L g over the practical range of tapes. For laminar 
flows the correlated pressure drops are weakly dependent on twist HID or HJD, yet for turbulent flows the pressure 
drop data are strongly dependent on twist. In most cases dpldz —> AP/L for single tapes in tubes. So in terms of 
Reynolds number the packed-bed Reynolds parameter .Y, - r! , un becomes 


Re 


2e 


Smithberg-Landis 


GqD p 


(l-£)|d 


= —Xt 


Ergun 


(A17) 


where 
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(A18) 




Ergun 


G 0 D p 
(1 — e)(x 


( f ^ ) Smithberg -Landis 

and the generalized relation can be expressed as 




v 

Ergun ^ Ergun 


(A19) 


y y 

^Ergun^ Ergun 


= 70.9 + 



(A83) 


From extrapolating the data g 2 (H/D, k ) -+ 0.014 and represents a lower bound of the data as illustrated in figure 3 
and labeled “single twisted tape.” To illustrate the dependency of the turbulent flow data on HJD, we normalized 
the Ergun friction factor by using equation (6). The modified Ergun friction factor is shown in figure 4 as 


(^Ergun) 


modified 


^Ergun 



f D ] 

2" 

1 + 




l Ho) 



(ii) 


and tends to follow the simplified form 


A+rgun ^Ergun 70.9 + 0.014X Ergun 


(A83a) 


TEST FACILITY 


The test facility consisted of an oil tunnel, flow system components, video equipment, a laser, lens systems, data 
recorders, and a test configuration (fig. 5). The test section consisted of 48 twisted (spiral) Lucite tapes assembled 
into a bundle and placed into a clear Lucite tube 2.54 cm ( 1 in.) in diameter. The tape bundle simulated, for 
example, a system of catalytic reactive surfaces, heat transfer augmentation surfaces, or flows through strata of 
porous media. The twisted tapes were made from clear, polished Lucite sheet 0.152 cm (0.06 in.) thick cut into strips 
0.318 cm (0.125 in.) wide by 17.8 cm (7 in.) long (nominal dimensions; measurements are discussed in appendix A). 
These strips were gripped 0.635 cm (1/4 in.) from each end and twisted with three complete twists (figs. 2(a), (d), 
and (f)). Figure 2(b) is a cross section of the actual packed bed, and figure 2(c) illustrates an ideal packed bed. A 
flow screen with 0.08-cm- (0.032-in.-) square mesh made of 0.023-cm- (0.009-in.-) diameter wire was placed across 
the tube inlet and attached to a square support (figs. 2(d) and (e)). For Borda inlet flows the screen restrained the 
axial movement of the twisted tapes and served as a screened orifice inlet when the test section was reversed. In the 
latter case, the twisted tapes were restrained 0.95 cm (0.375 in.) downstream by the visualization mirror. Some 
movement of the packed bed was noted. The movement would slightly disturb the inlet and exit flows but was not 
expected to alter the developed flow field. 

The assembly (fig. 2(d)) was then placed into a closed-cycle oil tunnel ( fig. 5) with a square support fabricated 
to retain the tube in the flow field and block the remaining tunnel cross section ( 152.4 cm by 7.62 cm by 7.62 cm; 60 
in. by 3 in. by 3 in.). Both the upper wall and the viewing port walls of the tunnel were Lucite (fig. 5). The index of 
refraction of the oil matched that of the Lucite, and magnesium oxide particles were used as flow tracers. The flow 
field was visualized by using the Full Flow Field Tracking (FFFT) method (Braun et al„ 1988). Laser light sheets 
illuminated two-dimensional sections of the tunnel along the flow path. Transverse visualization was accomplished 
by placing a mirror in the tunnel downstream of the 48-twisted-tape bundle and at approximately 45° to the flow 
axis. A second mirror placed above the tunnel (not shown in fig. 5) projected the view to the television camera. 

Cross tunnel traverses were accomplished by small rotations of the mirrors. These traverses provided insights into 
the three-dimensional nature of the flow field. 
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The coherent-beam, continuous-wave, argon-ion laser was directed by micrometric adjustable mirrors through 
two cylindrical lenses positioned at 90° to each other and through the Lucite tunnel window and into the test section 
The light sheet was approximately 0.01 cm (0.004 in.) thick, and the flow' was seeded with magnesium oxide flow' 
tracers. Micrometric adjustments controlled scanning of the light slices across the test section and provided a three- 
dimensional visualization of the flow field, which was videotaped at 30 Hz. 


EXPERIMENTAL RESULTS 
Visualization 

Figure 6 represents the axial flow field along the centerline of the packed bed of twisted tapes. Flows across the 
inlet were markedly influenced by the twist and packing of the bed. At the inlet the vena contracta normally found at 
the inlet of an open tube was limited to a minor region near the tube wall entrance (fig. 6(a)). Beyond this entrance 
region the flow was rapidly entrained into minor spiral perturbations (HJD 0 - 15.8) within less than 0.1D. How- 
ever, for flows near adjacent solid boundaries there is clear evidence of local spiral flows within the field of focus. 
Flow's in the boundary layer of each twisted tape tended to spiral as expected (fig. 6(b)). Thus, the flows tended to 
follow typical boundary layer flow patterns: no slip at the surface and a region of viscous flow closely aligned with 
the surface topology blending into a region of fully developed flow's with minor perturbation of the streamlines. As 
the pressure drop was increased, these regions became less distinct in that the boundary layers became thin and 
difficult to visualize. There still persisted a region perturbed by the spiral nature of the surface, but it became less 
distinct as the flow velocity increased. 

By aligning the laser sheet transverse to the flow we obtained a circular cross section of the flow. The packed 
bed of twisted tapes appeared as small, randomly oriented, rectangular blocks in the flow field (fig. 2(b)). There 
appeared regions where the flows w'ere jetting between adjacent blocks, regions where the flow's were predomi- 
nantly clockwise or counterclockwise, regions where the flow's were upward or downward, and wall boundary layer 
flows at both the “block” surfaces and the interfaces between the tube and the twisted tapes (fig. 7). Notably absent 
were local vortex flows, with the possible exception of one configuration. Although jetting should produce vortices, 
the three-dimensionality of the flow may have precluded observation. It appeared that engendered vorticity was 
captured within the axial spiral flow field, although not readily visualized. 


Packed Bed of Twisted Tapes 

The pressure drop and flow data parameters in table 1 were plotted in figure 8. The errors in the difference in 
static pressures upstream and downstream increased at lower Reynolds numbers, and those data should be consid- 
ered suspect. It is evident that the Borda inlet configuration with a screen had a slightly higher flow resistance than 
the orifice configuration with a screen and that the orifice configuration without a screen had the lowest resistance. 
Consistent sets of parallel lines could be passed through each set of data in figure 8. Shown for reference are the 
Ergun (1952) porous-media model and the Date (1974) single-twisted-tape model. In general, the data were low'er 
than the Ergun model over the range of data taken. The average fluid temperature was 22 ± 1 °C, the viscosity was 
0.91 poise, and the density was 0.9 g/em 3 . A suggested general form for N twisted tapes in a tube, following 
equation (A83), is 


X&gun pErgun - 45 + 0.009 X Ergun ( A85a) 

which is considerably below that suggested for flows in porous media 

%-gm^Ergun =150 + 1 .75X &gun (A84) 

suggesting significantly less flow resistance for the packed bed of twisted tapes. 

For flows through packed fibrous beds the resistance for parallel fibers is about one-third that for perpendicular 
fibers 


NASA/I’M — 2002-208914 


6 



(%gunX&gun) paraMel 36A- 3 ( y E r gun^E rg un) pel . pendlcular (A 

These expressions agree with those presented by Hersh and Walker (1980), but an average of the two forms more 
closely approximates the data of Sullivan (1941) 


y Y 
^Ergun^Ergun 


46.6 + 58.8 


= 52.7 


and is similar to the reported (laminar) data (table 1 and eq. (A85)). 


(A109) 


Single Twisted Tape 

Parameters V Er , nin and Ip from the selected inteipolated data of Date (1974), Koch (1958), Gambill and 
Bundy (1962), and Smithberg and Landis (1964) are plotted in figure 3 for the range 0.56 < HJD < <». The depen- 
dency on HJD is evident and is bound by the porous-flow loci (eq. (A84)) to that of a single twisted tape at high 
Reynolds numbers (eq. (A83)). This conclusion is based on the data of Koch (1958), which were taken from figure 7 
of Smithberg and Landis (1964). 

The results of Sparrow and Haji-Sheikh (1966) for laminar flows can be expressed as (see page 17) 

lErgun^Ergun = 70.94 (A39a) 

Also plotted in figure 3 are the data for 48 twisted tapes in a tube along with equation (A84). 

Figures 9 and 10 represent the loci of 48 twisted tapes in a tube, a single twisted tape in a tube, porous-media 
flows, and interpolated turbulent flow' data for a single twisted tape in a tube corrected for swirl velocity. 


VIDEOTAPE RECORDING 


The complexity of the flow field, whether virtual or experimental, became vivid through the videotape record- 
ing, which is included as .avi and .mov files on the supplement CD in the printed version of this report. Visualiza- 
tion of the flow field reveals flow threads, wakes, stagnation zones, and the influence of the twisted-tape interfaces. 
The flow threads can be observed during a scan of the flow field from the front to rear lateral walls. These flow 
threads weave through the packed array of twisted tapes in the bulk flow direction. Details of the flow boundary 
layer close to the wall and progressing through the packed bed to the opposite wall are revealed. The video can also 
be used to determine quantitative experimental information, such as the flow velocities, by using the FFFT tech- 
nique (Braun et al„ 1988). 
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SUMMARY OF RESULTS 


The flow experiment consisted of three principal elements: an oil tunnel 7.6 cm by 7.6 cm (3 in. by 3 in.) in 
cross section, a cylindrical tube containing the packed bed of twisted tapes in an arbitrary array, and a flow charac- 
terization methodology. Full Flow Field Tracking (FFFT). The indices of refraction of the oil and the test matrix of 
twisted tapes were closely matched, and the flow was seeded with magnesium oxide particles. Planar laser light: 
provided a two-dimensional projection of the flow field, and a traverse simulated a three-dimensional image of the 
entire flow field. Flows were observed near the inlet of the cylindrical tube housing the bundled array of twisted 
tapes, at the interface between the tube wall mid the twisted tapes, and within the bundle of twisted tapes. 

The flow field was three-dimensional and most complex to describe. The most: prominent finding was flow 
threads. The axial flow appeared to spiral along the twisted tapes within the confines of a virtual distorted cylindrical 
boundary. The flow field appeared to be simulated by a packed array of very thin virtual cylinders, with the 
exception of the spiral effect due to the twist. The effects of random packing and bed voids created vortices and 
disrupted the laminar flow but minimized the entrance effects of the unpacked tube. 

The results of several investigators for flows in geometries with a single twisted tape were analyzed. These 
results are related to the Ergun model in appendix A. (Symbols are defined in appendix B, and a comprehensive data 
table is given in table 1.) The single-twisted-tape results of Smithberg and Landis (1964) have been used to guide 
the analysis. The data for 48 twisted tapes in a tube were correlated by using the Ergun model for flows in porous 
media. The pressure drop and flow data for the three geometric configurations (Borda and orifice inlets with 
downstream restraining screen and orifice inlet without screen) have distinct flow characteristics differing up to 
13%. The averages for the combined data sets were lower than the Ergun model by a factor of 3 for the packed bed 
of 48 twisted tapes in a tube and by a factor of nearly 1 .6 for a single twisted tape in a tube. These results suggest a 
lower flow resistance for a packed bed of 48 twisted tapes in a tube than for either porous-media flows or single- 
twisted-tape flows. Further investigations including different geometric configurations and computational fluid 
dynamics analysis are suggested. 
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APPENDIX A 


RESULTS OF SEVERAL INVESTIGATORS SCALED TO ERGIJN MODEL 

Because we really do not know how to correlate the data between a single twisted tape and multiple twisted 
tapes in a cylinder or tube (figs. 1 and 2), we will begin with some familiar definitions of friction factor and 
Reynolds numbers and arrive at the form similar to that developed by Ergun for porous-media flows. 



fp u 2 L\ 

l 2 ) 

, d h , 


(general Fanning friction factor) 


(Al) 


R e=P^L 
P 

and for now, let p = Constant and CL = 4 /, as used, for example, by Sparrow and Haji-Sheikh ( 1966), Hong and 
Bergles (1976), and Gambill and Bundy (1962). 

Now let us define some packed-bed parameters. Let the average or superficial velocity within a packed bed of 
one or more twisted tapes be 



Uo 

£ 


where U is the velocity in the tube without tapes. The bed porosity is defined as 


(A2) 


£= 1 - 


Ntw 


D A 


(A3) 


where N is the number of twisted tapes, t their thickness, w their width, and D n the virtual twisted-tape diameter 
(fig. 2(c)). 



(A4) 


We can now' define the characteristic length of the packed bed in terms of sphere diameter. Note that D H ~ 
4AIS, where A is the cross section of the flow area and S the wetted perimeter. For a tube D H = D tube and for a 
uniform bed of spheres 6VJA S = 6/a v = D -» D here (Bird et al., 1960). This factor of 6:4 or 3/2 will become a 
scaling parameter for the Ergun relation. The characteristic length of the packed bed is 



(A5) 


where R h is the ratio of the bed cross section available for flow to the wetted perimeter, w'hich is equal to the ratio of 
the volume available for flow to the total wetted volume and is equal to the bed porosity divided by the ratio of the 
wetted surface to the bed volume. The specific surface area is 


a v = 


a 

1 — £ 



(A6) 
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where a v is equal to the ratio of the total sphere surface to its volume and is equal to the ratio of the wetted surface to 
the solid volume. 

Combining equations (A5) and (A6) gives the characteristic length of the packed bed as 



6(1 -e) 


(A7) 


In terms of one or more of N twisted tapes contained within a cylinder or tube of diameter /). and with equation 
(A3), 


‘ kD 2 


D h =- 


- Ntw 
> 


nD+2N(t + w) 


D£ 

n 2 N(t + w) 
kD 


(A8) 


When equation (A8) is substituted into equation (A7), D p becomes 


D, 


3D f1 t 

T (1 ~ e) 


1 + 


2 N(t + w) 
kD 


(A9) 


Note that in figure 2(c) the dark shaded areas A, B, and C are within the bounds defined by the twisted-tape 
width w and thickness t, the tri-circular loci with radii of DJ2, and the triangle with vertices 1,2,3. The area 
B + a + c = wt / 4 , the area C - c + b = wt / 4 , and the area A - a-b = 0. Therefore, the dark shaded area (solid area 
of the tape) within the triangle is will. Continuing with triangle 4-2- 1 , C t - Cj + bj = wt / 4 , Bj + Cj + a = wt / 4 , 
and A t - tq - a = 0 . Similarly, for triangle 2-5-3, C 2 - c + b 2 = wt / 4, B 2 + c + a 2 = wt / 4, and A 2 - b 2 - a 2 = 0. 
Summing these areas gives the total solid twisted-tape area within the confines of the hexagon. Extending these 
results provides a generalized form for N twisted tapes within the confines of a hexagonal space that approximates 
that of a circumscribed cylinder or tube. This generalized form is an ideal model for packed beds that is expressed 
only in terms of tape width and thickness. 


i- 2 A— 9 

£ model = T 2 ( A1 °) 

In Ain V3 C +w z 

2 0 2 0 

and for w = 2 1 (see also the section Some Sample Calculations at the end of this appendix) 

e= 1 = 0.538 

5V3 

Comparing the model to the experiment where 48 twisted tapes with average width and thickness <(w) = 0.3234 cm 
(0.1273 in.) and (t) = 0.1275 cm (0.0502 in.) were contained in a 2.54-cm- (1.0-in.-) diameter tube gives an average 
porosity, from equation (A3), 


/ \ 3 192(0.3234)(0.1275) 

VW- 1 27t tt(2.54) 2 


(All) 
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which indicates that the experimental test tube was well packed. For a single twisted tape in a tube (Smithberg and 
Landis, 1964), where t = 0.0559 cm (0.022 in.), w D g = 3.51 cm (1.382 in.). 


e Smithberg-Landis - 1 — - 0.9797 (A12) 

KL > 0 

A problem arises because this formulation is independent of twist HID. From structural mechanics the twisting 
of tapes is considered to follow linear mechanics with little change in length, unless plastic deformation with 
subsequent annealing takes place. As a result of this assumption the bed porosity is not a function of H. However, 
the flow path is a function of twist HID, and consequently, the characteristic length must be modified to reflect H. 
The spiral path is defined in terms of the parameter (p as 

jc=rcos(p; y = rsintp; Z=/ttp (A 12a) 


and the path-length amplification ratio becomes, for the number of 2it twists n v 


L _ 2nn,r 
7 7 


1 + 


2mi t r ) 


(A 12b) 


For very large HID (slow spiral), n, 1 and L/Z 0 -> 1. For very small HID (rapid spiral), n t -> M » 1 and UZ 0 -> 
2j iMr/Zg. This amplification is not strong enough to account for the pressure drop increases due to the twist. 


Models for Single-Twisted-Tape-in-Tube Analogy 


Smithbera and Landis (1.964) . — Smithberg and Landis (1964) used the local average velocity (superficial bed 
velocity) in their correlations. Consequently, their friction factor and Reynolds number relations can be scaled 
directly in terms of the packed-bed parameters T Er „ im and X El . gU[1 . Substituting equations (A2), (A7), and (A1 1) into 
equation (Al) gives the pressure drop parameter L Ergun in terms of the Darcy friction factor as 


4/ 4/smithberg-Landis 


[ A/ ’ 1 


1 

'pA P' 

e 3 ) 

( d \ 

U P 

^0.5pn 2 , 

U oj 

“ 3 

I G 0 2 J 

l 1-e J 

v L 0 j 


(A13) 


./Simthberg-Landis ^Ergun 


where 


G 0 =pU 0 


(A14) 


(A15) 


' i ■[:' un 


'pAP' 

{ £3 1 

1) 

p 

\g 2 0 ) 

i 1_e J 

K Lo j 


(A16) 


Here L 0 represents the straight-line distance between pressure laps. The twisted-tape length L varies little from L 0 
over the practical range of tapes, and for laminar flows the correlated pressure drops are weakly dependent on twist 
HID. However, for turbulent flows the pressure drop data are strongly twist dependent. In most cases dpldz -» AP/L 
for single tapes in tubes. Therefore, in terms of Reynolds number the packed-bed Reynolds parameter X Ergun 
becomes 
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Re 


2e 


Smithberg- Landis 


G„D„ 


(l-e)p 




Ergun 


(A17) 


where 


X, 


Ergun 


GqD p 

(l-e)|i 


(A18) 


2e 

(/ Re )smit h berg-Landis = ~ ^Ergun^Ergun 


(A19) 


For turbulent flows the F Ergun X Er product is dependent on both Re and HID. Smithberg and Landis (1964) 
provide a simple expression for the Fanning friction factor. It is important to recall that herein we used HID - 180° 
twist, whereas Smithberg and Landis used HJD = 360° twist, where H and H a are measured along the axis parallel 
to the tube centerline, 2 


H, 


O 

D 


= (-) - 2^1 

V D /Smithberg-Landis \D J 


(A20) 


/s niithberg- Landis 


0.046 + 2, 


{ D 


- 1 . 2 ' 


Re - ” 


(A21) 


where 


n = 0.2 




Adjusting the constant 2.1 and the exponent -1.2 slightly gives a better fit to the data, and the normalized form 
becomes 


' (H 0 \ 


£ 

D J_ 


yfo j 


(*-r 


! 1 + 1051 

^ Smithberg- Landis 

where/ = 0.046Re a2 . In terms of equation (A21), equation (A19) becomes 


Re 


-0.34 /JH 0 /D 
Smithberg -Landis 


(A22) 


y y = 

1 Ergun ^Ergun ^ 


An*: 

8 °\ D 


(. fo Re ) 


0.207 

Smithberg- Landis g 


So 


Re 


0.8 

Smithberg-Landis 


(A23) 


After substituting equation (A17), equation (A23) becomes 


Y \ irgun A[ j-gun 0.158 


0.8 


( H Y 

A/J 


Aj 


0.8 

Ergun 


(A24) 


“Yet the number of 180° twists will be twice the number of 360° twists for a fixed tube length (i.e., 80 ° twists = 2w 360 o twists ), 

and some authors use twist count rather than measured values of H and H g 
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where 


- r *y 


[i + iosfM_i 

x— 1.15 

( Xi: \ 

o £ Ergun 

A l)J 


V D 

J 

l 3 J 


At very high Reynolds numbers surface roughness will promote eddy bursts near the walls and separation effects as 
the flow attempts to track the twisted tape, implying that (figs. 8 and 10) 

y Ergun — > Constant — » 0.014 for e — > 1 and Re » 2000 (A25) 

where 0.014 is the coimnercial rough-tube equivalent (see eqs. (A54) and (A55)). For two-phase flows excessive 
vapor generation, or holdup, would be anticipated. 

Hong and Bergles (1976) . — A similar set of scaling parameters can be developed for the data of Hong and 
Bergles (1976): 


Re 


Hong- Bergles 


puD 


r G 0 D p \\ 


D 


From equation (A8) with N = 1 and (w + /) — > A 



f , . t + w 3 


A 23 

D 2 

1 + 2 


1 + — 

K D 

2 

— ^ — 
3 

K 

D p ~ 3 

1-8 

V J 

1-8 
v y 


(A26) 


(A27) 


Re 


Hong-Bergles 




1,09 

e 


A, 


Ergun 


(A28) 


C 


/, Hong-Bergles 


AP 

f— 1 = 2 

'pAP' 

D 

p 


^0.5pr< 2 , 

U" J 

{ L J 

KJ 


e = 23 Ergun T 



fl + 2 f + W ] 


fl + -l 

K 

2 

nd 

4 

3 

8 

— ^ — 
3 

8 


V ) 


V. y 


'Ergun 


(A29) 


where LJL ---> 1 . With t - 0.046 cm, D a = 1.02 cm, and w — » D g . 


0.046 _ nn ._ 

£ Hong-Bergles ^ — v.yQj 

1.02 — 

4 

Equation (A10) would give 0.945, so w — » D a is a good approximation. 

O/.Hong-Bergles — 2.3141g r g Un 


For laminar flows in a half-tube configuration, equivalent to a tube with a single twisted tape. 


(A30) 


(A31) 
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c /Re ) 


_ 8 

Hong-Bergles 9 


K 


^2 


^Ergun^Ergun 183.6 


(A32) 


^Ergun^Ergun - 68.6 


We will later show the relation to the work of Sparrow and Haji-Sheikh (1966) as 

2 

(Cr Re) = 183.6 = f^-) (CrRe) 

V 1 /Hong-Bergles \ D ) ' 1 'S 


/Sparrow-Haji-Sheikh 


(A33) 


(A34) 


where D H !D = 0.5682 is found from equation (A8) by using the geometry of Hong and Bergles (1976), where t = 
0.046 cm, w = 0.97 cm -> D 0 , D = 1 .02 cm, and (D H ID) 2 0.334. 

Date (1974) . — For the work of Date (1974) the scaling factors become, for (t + w) -» D, 


Re 


Date 


pud 


r P UpDp 


D 

K 

£ 

V J 


K_ 

£ 

V ) 


X, 


Ergun 


(A35) 


/Date 


J _(dp_ 

2 Vdz 


D 

2 

P M 


\_( dp 

'l{ ~dz 


D„ 


v „ 3 


l p4 , 


D 

K D PJ 


*4 


'i + i' 

n 


'Ergun 


(A36) 


(/Re) 


_ 2 
Date — g 


1 + ^ 
n 


£ 

\ / 


l/rgun^Ergun 


(A37) 


For laminar flow in the half-tube configuration, or for a tube with a single twisted tape, and e-»l (see Weigand, 
1948, and eq. (A76)) 


/Date Re D ate= 42.19 (A38) 

^Ergun^Ergun = 20.9 (A39) 

which is in good agreement with equation (A33). For e = 0.934, F Erguff X Ergun =61.8; and for e = 0.9797, 

^ Ergun^Ergun ~ 68.0. 

For Re Date < 150, tape twist HID has little influence on the Fanning friction factor. For 150 < Re Date < 1000, 


/Date /bate-half-tube A/bate 
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(/ Re )p a te 


/Date-half -tube Re Pate 


1 + - 


A/Pate R( -Pate 


/bate- half -tube Rl -Pate 


( 

.19 1 
V 


= 42.19 1 + 


Afoate Re Pate 
42.19 


1 + A/bate 

f , 

2 1 + _ 

2 K 

^Ergun 

2 

fl + -l 

K 

3 8 

42.19 

9 

8 


V ) 



V J 


V Y 
' Hr gun / M’igun 


(A40) 


and for e -+ 1 


^Ergun A rgun 70.9(l + 0.026A/ Date X Ergun j 


(A41) 


and for 150 < Re Date < 1000 


A/bate = 0.088 


-1 

H /Pate 



(A41a) 


For Re Date > 1000, the Srnithberg and Landis ( 1964) expression, or modified expression, for the Fanning friction 
factor can be converted to/ Date . Recall equation (A2 1): 


/s mithberg- Landis 


0.046 + 2.1 — — — 0.5 


where 


n = 0.2 


1 + 1.7 


(Ho. 
V D 


ivr 


-1.2 


Re - ' 1 


and the conversion for Re Date > 1000 follows as 

(/ Re )pate 


( D A 2 
V °H ) 


0.046 + 2.1 


"2 H) 
~D) 


n-1.2 


-0.5 


>Re 


1 -m 


Date 


where 


(A21) 


(A42) 



r -1 

-0.51 


(H\ 


+ 

2 — 

V D /£) a t e 

1 


The DID h = 1.66 given by Date (1974) should be nearly 1 .6366. Date gives no values for tape width or thickness, 
and estimates of £ will give a range 61 < (L Emun A|; rgun ) Date < 7 1 for laminar flow in a half-tube or in a tube with a 
single twisted tape. Also note that Date's calculations underpredict turbulent friction factor data by 30%. 

Gambill and Bundy (1962) and Gambill et al. (1961) evaluate swirl-flow heat, transfer along with isothermal 
data on friction coefficients. The relation given 
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where 





(A43a) 



(A43b) 




izD-2t + 2D 


(A43c) 


and where y is the number of tube diameters per 180° of twist (2 nfl), D e the equivalent or hydraulic diameter (over 
the practical range 5/4 < y < », D g ■■■■■> 8 the surface roughness, f s the Darcy swirl friction factor, and/ 0 the 

Fanning friction factor. The reduction of the analysis parameters to Ergun parameters takes the same form as that of 
Sparrow and Haji-Sheikh (1966): 


fa~Cf ~ p| 


-dp\2 
dz ) 3 




where — — » Constant > 0 


dz 


(A44) 



However, the data presented in figure 2 of Gambill and Bundy (1962) follow the parameters of Hong and Bergles 
(1976): 



(A28a) 



(A29a) 


& C /,Hong-Bergles 2314} Ergun (A3 la) 

Equation (A43) represents the best fit of a large data set from several investigations, yet at Re g = 10 000 and 8 ID = 
ox 1 () 4 the ratio of friction factor data at y =1.12 to those at y = 0.28 (an unusually tight twist) is nearly 17. From 
equation (A43) the ratio is 6. Although this difference has not been resolved, it is important to recognize that at a 
given Reynolds number and value of y (noting that HID = y ID) the effect of surface roughness is very important 
because both increase the friction factor. 

Sparrow and Haji-Sheikh (1966) . — For the work of Sparrow and Haji-Sheikh (1966) on flow and heat transfer 
in an arbitrarily shaped tube, the scaling factors for the half-tube configuration become 
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= 63.06 


(A46) 



After substituting equations (A 2) and (A7) and recalling our assertion that dp/dz —> AP/L, equation (Al) becomes 


C f 



2 

3 



4 

3 


' 1 1 'Min 


where 


-dp 

dz 


* 


Constant > 0 


(A44a) 



^Ergun^Ergun — 70.94 


Sparrow and Haji-Sheikh (1966) also estimate entrance losses as 


Ci 


r _z ' 

\°HJ 


+ K = ■ 


A P 

0-5p«ba 


(A45a) 

(A39a) 


(A47) 


where 



u 

^bar > 


dA 


and for the half-tube configuration 


K 1 .463 (A48) 

The agreement between Sparrow and Haji-Sheikh (1966) and Date (1974) for the half-tube or single- twisted-tape 
configuration is not surprising as they arc within \.% of the values of Weigand (1948), which we discuss in the 
section Torsion-laminar flow analogy . 

Bird et al. (1960) . — For laminar flow in a tube (limit e — > 1 ), Bird et al. (1960) give the average flow velocity as 


u = ■ 


APRl 

2\iL 


Substituting equations (A2) and (A7) and multiplying both sides by p “U give 


Cl = 0.51 


pAPe 3 

d 2 \ 

u v\ 

M 


36(1 -e) 2 

y Lq J 

{ L ) 


(A49) 


(A50) 
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1 = 0.5 


Y Y 
^Ergun^Ergun 


(A51) 


36 


V } 


^Ergun A Ergun = 72 T L where ~T~ 1 (A52) 

which is in good agreement with the results of Date (1974), Sparrow and Haji-Sheikh (1966), Hong and Bergles 
(1976), and Weigand (1948) for laminar flows in the half-tube or single-t wisted-tape configuration. 

From data for flows in packed beds Bird et al. (1960) give 2 IJL 0 = 25/6, resulting in the Blake-Kozeny equa- 
tion. Kaviany (1995) gives 2 UL 0 = 5. From these results we anticipate that the bed correction factor for laminar 
flows in a configuration with N twisted tapes in a tube will be largely independent of the twist for 



'l + l' 

K 


X, 


Ergun 


<150 


V J 


where X Ergun < 400 / 3 and e -» 1 


(A53) 


with some correction for flow path length 2 UL 0 . However, the correction factor will not be as strong as that for a 
randomly packed bed, where < 10 is considered laminar. 

At high Reynolds numbers the results of Smithberg and Landis (1964) show significant effects of twist HJD as 
discussed earlier (eq. (A21)). At even higher Reynolds numbers the surface roughness will promote both eddy bursts 
near the walls and separation effects. For commercially rough tubes (Schlichting, 1955) and in the limit: e-»l (i.e., 
without a twisted tape in the tube) 


( 10 ^ Re Wuradse(1933) = 400/( f ) X Ergun «^FErgun f ^Ergun = l- 27 (A54) 

^Ergun Constant -¥ 0.014 for e -> 1 (A55) 

and 0.014 is the commercial rough-tube equivalent. From the data of Smithberg and Landis (1964) the trends for the 
single-twisted-tape or half-tube configuration, using D H , appear to approach equation (A55). 

Torsion-laminar flow analogy . — Weigand (1948) provides solutions to the problem of torsion in prismatic 
members where the function J{x,y) satisfies (see also tig. 2(g)) 

V 2 / = - 1 (A56) 

with f b = 0 on the boundary of the cross section. This function gives the torsion constant J d for a member of cross 
section A, where 


Jd 



(A57) 


The analogy between fully developed flow in a tube and torsional problems starts with the Navier-Stokes 
equations 


F-Vp / p + vV 2 u 


(A58) 
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assuming steady flow where inertia and body forces are small and V/; = Constant. It then follows that 


V 2 u = -1 


with r.< bar = 0 on the boundary where 


and 



(A59) 


(A60) 


(A 61) 


The equations describing this type of flow are in a form similar to those for torsion in a prismatic bar as 
described by Weigand (1948) and pointed out by K.C. Cheng in a discussion to the work of Sparrow and Haji- 
Sheikh (1966). This leads to 


— V 2 ««-1 = V 2 / 

4C„ 


(A62) 


/ = 



u = u 

' -dp N 
dz 

l 4 J 





l y 


(A 63) 


M bar — 


JJ 


udA 


-\U 


-dp i 


' -dp ' 

dz 

dA J d 

dz 

, u 

~A~~ 

M- 

k j 


c ) 


\a~ 


(A64) 


Now the coefficient of friction and Reynolds number are related as 


C/ Re = 


-dp ' 

dz , 


ft 


0.5p«bai 


pMbar ft 


dz ) 
0.5Ubarlt 


(A65) 


Substituting for « bar for fully developed laminar flows, where D H = D e , gives 


c /Re = 


8 Aft- 


’d 


(A66) 


Substituting for 
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4A 


D„ = - 


Perimeter 


gives 


where 


C/Re = 


128A 


J d (Perimeter)" 


A = kR a 


1- 


f 2a -sin 2a ^ 
l 2jt )_ 


(A67) 


(A68) 


(A69) 


Perimeter = 2kR\ 


1 - 


a - sin a 


(A70) 


For circular tubes (x = Jt - (a = 0) -> n) Weigand (1948) gives in his table IV 


K = -4- = 1.571 = — 
R 4 2 


(A71) 


and 


C ; Re = 64 


(A 72) 


For semicircular tubes [x = Jt - (a = «/2) -4 rc/2J Weigand (1948) gives in his table IV 

k = —4- = 0.298 
R 4 


(A 7 la) 


h 


CV Re = 16tr 0 

4 / O \Z 


/ 9 V 

i+- 


= 62.97 


(A 73) 


(C /R e), 


Sparrow -Haji- Sheikh (1966) 


= 63.06 


(A46a) 


Although the (.(Re values of Sparrow and Haji-Sheikh (1966) are monotone with a, those of Weigand (1948) are 
not. However, all are well represented by 


C f Re = 63 = 


128 A 3 

/^(Perimeter) 2 


(A74) 


and would be a good engineering approximation for laminar flows in tubes of various cross sections (e.g., semicircu- 
lar). For Date's work, using the preceding relations. 
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Cf Re = 4f\ 


= 4 f**(T = 63 


and for the twisted tape under the following conditions: 


H k D , 2 

y = a = — ; — = 1 + — 

I) 2 l\ it 


/-Re,- =42.19 


Generalization . — For tubes with single twisted tapes 


^Ergun^Ergun 70.9 + 


Wf )]0 


where from Smithberg and Landis (1964), for X &gu[) > 2600/3, 


—11 = 0.046 + 2./— -02 

D)\ \ D 


n = 0.2 1 + 1 


-(f)' 


and for 400 < 3X Freun < 2600 


= 0.088 


— + 0. if — 1° 1 1 log 10 
H \H ) 610 3 e 150 


Some consolidation in the turbulent flow regime follows from Smithberg and Landis (1964): 


-^— = 0.9 + -^ 

fH/D^oo f±o_ 


. _ 0.046 

J HI D ^> oo - Re0 2 
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. Thus, 


and is thought of as the nonideal flow path length L for a particle to travel a bed or tube length of (7 7; ) ///;) l 
g^H/D) can be modified as 



(A80) 


Lopina and Bergles (1969) . — Of equal interest is the simple relation given by Lopina and Bergles (1969) for 
turbulent flows in smooth tubes (+20%) for water and air for 5xl0 3 < Re < 4.5xl0 5 and < H/D < ^/85 : 



- 0.406 


(A81) 


where f 0 is the isothermal friction factor for an empty tube,/, = 0. 046/Re 0 ' 2 . Reducing the results of Lopina and 
Bergles (1969), with equation (A9), gives 


/Re = 


( ARAP \ 

( GD, \ 

{ D p e l-e 

l 2 LG 2 J 

{ H J 

[Dp e l-e y 


£ 

e 


(A81a) 


/ 0 Re 


' 1 ? 
\fo J 


— ^Ergun ^Eiguii 0-6 1| 


(!) 


- 0.406 


1 + - 


(A81b) 


where 


/ 


1 +^ 

K_ 

3e ? 


'Ergun 


V ) 


(A81c) 


Re = 



(A81d) 


The resulting fit to the data is lower than the data of Smilhberg and Landis (1964) by about 10% for HID = 2.48. So 
the agreement can be good and the simplicity may be of merit for a limited Reynolds number range. 

We also note that for large values of X &oun the data will become independent of X[; r>un as equivalent surface 
roughness (k s s 8 ID) plays a role. Thus, for large .Y,,- rL)[||| for a given geometry 


^ Ergun 



= Constant 


(A82) 


or in generalized terms 


y y 

^ Ergun J Ergun 


= 70.9 + 



(A83) 
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and from extrapolating the data 


^Ergun^Ergun — 70.9 + 0.014A Er g Un 


(A83a) 


Porous Media Flows and N Twisted Tapes 

For porous-media flows Ergun (1952) combined the turbulent and laminar data, noting that at high Reynolds 
numbers the friction factor approached a constant as is consistent with flow in a fully roughened tube (k s > 0). The 
result was a linear sum of the Kozeny and Burke-Plummer equations, and as illustrated in Bird et al. (1960) the 
generalized form may be written 


7"t r-Min ^1 i -'iin = 150 + 1.75X &gun (A84) 

And if the Ergun data were for a single insert, and following the development on page 199 of Bird et al. (1960), 


P 


2 V 2 

v 0 


= pAP 


f D \ 

f " 3 ] 


" p v 0 d p " 

U J 

I 1 " 6 ] 

l 72 J 



^Ergun ^Ergun 7 2 


which nearly coincides with the 70.9 of Lopina and Bergles (1969) in equation (A83), indicating that friction factors 
for flows in tubes with multiple twisted tapes should be much higher than measurements are showing. 

For tubes with N twisted tapes the laminar data are below, yet parallel to, the 7f Er£ , im T Er0un relation for flows in 
porous media. From the data in table 1 (which may yet require Prandtl number or viscosity corrections). 


-^Ergun ^Ergun ^ Onstant 45 


(A85) 


with a “suggested general form" for N twisted tapes in a tube following equation (A83), w'hich appears as 


•^Ergun ^Ergun 45 + 0.009 7f Ergun 


(A85a) 


From equations (A93) and (A95) for parallel flow's in cylindrical fibrous materials, with e = 0.61 for a 48-twisted- 
tape cylindrical bundle. 


^Ergun^Ergun ~ 36fc - 46.6 (A85b) 

suggesting that pressure losses in the twisted-tape bundles are less than in other porous-media flows. 

The combination of implied losses and those of, for example, cotton fibers, suggests a flow vortex structure that 
impedes the passage of fluids more than hair or glass fibers do. This structure does not suggest a direct relation to 
either heat or mass transfer. These implications remain to be investigated. 


Fibrous Buik Materials 

The work of Fowler and Hen el (1940) for flows through wads (e.g., wool, glass wool, cotton, rayon, kapok) 
provides a solution: 



(A86) 
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where 0.18 < k < 0.2, j a P = p, and the specific surface area a = x/g equals the ratio of the element volume x to the 
element surface c. Equation (A86) can be rearranged to the Ergun form: 



( e 3 ] 


(-dp" 

1 (l-e)p ,/ x ) 2 

( \ 

P 

dx 

A|irgun G 0 D p \C ) 

[i-ej 

Wo) 

D p 

v / 



(A87) 


where 


and 


For k ~ 0.2 


*Ergun 


" pAP 

r e 3 ] 

/) \ 

Go J 

i 1_e J 

K L O J 



^Ergun^Ergun — ISO 

as also recommended by Kaviany (1995). For k - 0. 18 as recommended by Fowler and Hertel ( 1 940) 


(A16) 


(A88) 


(A89) 


^Ergun^Ergun — ^00 (A 90) 

The work of Sullivan (1.942) for parallel fibers illustrates a distinct departure from the relation F Ergun X Ergun = 
Constant. Sullivan’s relations are similar to those of Fowler and Hertel (1940) and follow the same reduction to the 
Ergun form: 


D„ 



(A91) 


1 (l-e)p 

( \ 
5 

( e 3 ) 

Al'rgun G 0 D p 

k V 2 /) 2 

y K o°o LJ p y 

l 1_e J 



^Ergun^OErgun 36 k 


where the units of Sullivan (1942) are in the cgs system (dyne, g, cm, s). 


(A92) 

(A93) 


k = ^- (A94) 

S 

and p = (sin 2 <j)) avg , where <|) is the angle between the interface normal and the microscopic flow. For flows with 
parallel “cylindrical” fibers and »1 
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kq — > k L 


e 

3 


+ 0.8 


(A95) 


O 


1.02- e 


For flows with parallel cotton fibers and e < 0.85 


kq — > k Q = 2.5 


(A96) 


For £ > 0.85, k‘Q—>k 0 tends to follow k Q para u e i. 

The ^Eroun^&aun P r °duct for parallel flows in cotton fibers with e < 0.85 is about half that cited by Fowler and 
Hertel (1940) for packed wads (90 versus 180 or 200 depending on the value of k used for packed wads of fibers) 
and would agree with Fowler and Hertel for packed-wad flows where e -+ 0.95. The k used by Sullivan is the 
inverse of that used by Fowler and Hertel. 

Hersh and Walker (1980) revisited the work of Sullivan (1942) and provided a correlation over the range of data 
as follows, where the units are in the cgs system (dyne, g, cm, s): 

For flow parallel to fibers (0.1 < e < 0.985) 


2/ Re = 


APct 

liLU 0 


15.74(1 -e) 


1.413 


1 + 


27(1 -e) 3 =, 


'parallel 


(e) 


For flow perpendicular to fibers (0.7 < e < 0.992) 


2/ Re = 


APct 

iiLU 0 


= 64 


(1-e) 


3/2 


1 + 14.75(1 £) ^perpendicular ( e ) 


(A97a) 


(A97b) 


where the latter expression was developed by Davies (1952). 

After refitting the data of Sullivan (1942), the Hersh and Walker (1980) form may be re-expressed as 


' APd 2 ' 
v M -LUg ^ 


parallel 


E 

3_ 

1.02-e 


■ + 0.8 


n-ef - 1 

'A Pd 2 ' 

V e J 3 

. V-LUo , 


perpendicular 


(A98) 


And the Ergun form follows: 


(^Ergun^Ergun) parallel 36 k ^ (FErgun^Ergun ) perpendicular (A "' 

Although these expressions generally agree with those presented by Hersh and Walker (1980), an average of the two 
forms more closely approximates Sullivan’s (1941) data (see figs. 2 and 3 in Hersh and Walker, 1980). These 
modified forms (eqs. ( A97) and ( A98)) are presented herein as figures 1 1 and 12. 

Note that for parallel flows through cotton fibers the shape factor is 2.5, and for flows though fibers aligned 
normal to the flow (Davies, 1952) it is 3. The implication is that flows though fibrous materials (e.g., cotton) aligned 
with the flow do not differ significantly from flows through fibrous materials aligned normal to the flow. However, 
flows through fibrous materials still have three times the pressure drop as flows along parallel cylindrical materials 
when other parameters remain fixed. Sullivan (1941) found that the pressure drop through packed beds with fibers 
oriented perpendicular to the flow was double that of flows through parallel fibers. 
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Reducing the Hersh- Walker (1980) form to the Ergun (1952) form, noting that Sullivan gives (d) as follows 
from equations (A6) and (A7), results in 


£H n D P E 

4 h 6(1 -e) 


\ S o J 


1-e 


where 


(d)~ = 


4 2 D p 


S„ 3 


Substituting, equation (A97) becomes 


f \ 

( 3 \ 

_p_ 

£ 

Wo) 

1 1-e J 


2D, 


AP- 


■ = F{e)[iU c 


W; 


( „3 \ 


v !-eJ 


•^Ergun^Ergun 




For flow parallel to the fibers (0. 1 < e < 0.985): 


(A100) 


(A101) 


(A102) 


(A103) 


JErgun^Ergun = 35.415{[(l-e)-°- 587 e 2 ][l + 27(l-e) 3 ]} (A104) 

and, for reference, as e -> 0.27 the quantity in braces approaches 1. Equation (A104) can be represented by 
0.75 <e<l. 


(•^Ergun^Ergun) —1.6 0.71Og e (l e) (A104a) 

For flow perpendicular to the fibers (0.6 < e < 0.992): 

%gun^Erg„n = 144{[(l-e)-°- 5 ]e 2 [l + 14.75(1 -e) 3 ]} (A105) 

Equation (A 105) can be represented by 0.75 < e < 1. 

(^Ergmi^Ergun)^ 4 = 2.5 - 0.81Og e (l - e) (A 105a) 

and for this case the quantity in braces approaching 1 as e — » 0.5 is beyond the region of validity of the relation. 

From figures 1 1 and 12 (replotted figs. 2 and 3 from Hersh and Walker, 1980), the relations appeal - to be in 
good agreement, yet the sensitivity to small changes in e becomes paramount. Compared with the twisted-tape data 
(see table 1) for e = 0.522 and A = 48 twisted tapes in a cylindrical (parallel) bundle, equation (A104) gives 
^Eroun^Eroun = 58.8, which is nearly 7/3 greater than the data (see eq. (A85)). As the bed porosity approaches that of 
a filter or a particulate separator, the 2f Er „ un F Ergun product is not a constant but depends on the filter porosity, as 
noted by Hersh and Walker (1980) and Sullivan (1942). 
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Silverman and First ( 1952) reported data for edge filtration of 0.5-p.m dust and a baby-oil smoke generator with 
average 0.6-p.m particulates. The filtration materials were a variety of felts, flocked papers, and fiberglass media. 
They noted that rockwool, silica fiber, or aluminum oxide fiber could be used at high temperatures. For their 
correlation 


JErgun^Ergun = ^29 £ 3 (1- e) L4 (A106) 

Davies (1952) provides a theoretical foundation for several types of particulate separator, such as settling chambers, 
conical elutriators, inertial and cyclone separators, jet impingement mechanisms, precipitators, thermal separation 
mechanisms, and filters. Data for flow through fibrous materials follow: 

5fegun^Er g u„ =^(70)e(l-e) 1 ' 5 [l + 52(l-£) 1 - 5 ] (A107) 

The Davies (1952) results range up to three times higher than those of Silverman and First (1952) to several times 
higher than those of Ilersh and Walker (1980). 

Using the data in table 1 for 48 twisted tapes in a cylindrical (parallel) bundle, where e = 0.522, gives, 

^Ergun%gun = 58.8 from eq. (A104) (A108a) 

-^Ergun ^Ergun = 46.6 from eq. (A85b) (A108b) 

For the 48-twisted-tape data prediction 

XErgunlErgun = ^ \ ^ = 52.7 (A109) 

or less than 1.2 times that of the reported data (table 1 and eqs. (A85)). 


Packed Beds of Spheres 

Wentz and Thodos (1963) measured the pressure drop across packed (cubic, body centered, and face centered) 
and distended bands of five layers of 3. 1 2-cm- ( 1 .23-in.-) diameter spheres held in place by short wires in drilled 
holes and epoxy. Spheres in the distended models were separated to simulate bed swelling. Both sets of data were 
correlated by 


^Ergun = v ,ms 396 , I 2 550 < * Ergun < 64 900 (A1 10) 

TtErgun “ U2 

which is less than half that of the Ergun equation (A84) but over a larger range of X Ergun . As the bed length was 
short (five spheres), X El .„ im dependence may be attributed to as-yet-undeveloped turbulent flow. A problem combin- 
ing or extending this relation occurs with the laminar regime because a singularity occurs at X Eroun = 38.34. 

However, 


*Ergun%gun = 150 + 4.5X^ 8 g 4 un (A1 1 1) 

is a form that includes this high-Aj. un turbulent regime, yet is higher than the Ergun equation (A84) for the 
transition region 15 < < 300 and as much as 14% higher for the region 40 < X Ergun < 80 (fig. 13). The 
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turbulent-flow pressure drops through these packed beds of spheres have similar dependence yet are 30 to 40 
times larger than for a single twisted tape in a tube over similar ranges in The expression 

X&gunJfegun = 003 [ 150 + 4 -5X&g 4 un] (A1 12) 

provides a reasonable fit to the data of Smithberg and Landis (1964) and Koch (1958) in figure 14. 


Some Sample Calculations 

The basic parameters for 48 twisted tapes in a 2.54-cm- ( 1 .0-in.-) diameter tube are as follows: 

Fluid: 

p = 0.9 g/cm 3 
p = 0.9 g/cm-s 

Tapes: (w) and (t) represent average values of twisted-tape width and twisted-tape thickness measured at each end of 
the 48 twisted tapes (table 2). 

<w> = w = 0.3234 cm (0.1273 in.) 

(t) = t = 0.1275 cm (0.0502 in.) 

D 0 = (m- 2 + t 2 ) 0 - 5 = 0.348 cm (0.137 in.) 

Three full 360° twists in 16.5-cm (6.5-in), length provides an average twist ratio of 


16.5 

-^- = — — — = 15.8 
D 0 0.348 



L 0 = 1 6.5 cm 

Other parameters: (e exp ) = 0.61 is the average area-weighted porosity (0.593 < e exp < 0.628). 

£ exp = 0.61 


2 Z) (l-e) 

Dp ~ 2N(t + w) 


(1.5)(2.54)(0.39) 


1 + - 


n D 


1 + (2)(48) 


' 0.1275 + 0.3234 A 


= 0.231 cm 


(2.54)7t 


A n = 2.54' 


if n' 


= 5.067 cm 
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dV_ dV 09 

G n = — — = dt 60 _ 2 95 x io “ 3 y where V is in cm 3 / min 
° A 0 5.067 

AP = dP(6.8947xl0 4 ) ing/cm-s 2 , where P is in psia 


oDp 0.231 u ,, . . i 2 

— f— = G„ = 0.659G„ where G n is in g / cm - s 

~ \ . is / r\ o r\\ / r\ r\ \ is u a 


'*** ~ “ W 1(0.39)(0.9)J “ WUC 

v _ ( pAP )f Up Y e 3 \ (0.9)(0.61) 3 0.231 A P 


ErgUn [gZ \{L 0 ){l-e\ 0.39 16.5 G 2 


= 73.34x10“ 


locale = Y^~ + 1-75 
-^Ergun 



Geometry 

dV/dt, 

cm 3 /min 

A P, 
psi 

O, 

g/em'-s 

dP , 

g/cm-s 2 

^Ergun 

^Ergun 

^calc ^cal</ 

^Ergun 

Borda + 
screen 

5900 

2.54 

17.46 

17.5X10 4 

11.5 

4.22 

14.8 3.51 

Orifice + 
screen 

4754 

1.88 

14.07 

13.0X10 4 

9.27 

4.81 

17.9 

3.73 

Orifice - 
screen 

6767 

2.47 

20.0 

17.0X10 4 

13.2 

3.12 

13.1 

4.21 

Borda + 
screen 

6777 

2.86 

20.0 

19.7X10 4 

13.22 

3.6 

13.1 

3.64 

Orifice + 
screen 

6556 

2.48 

19.4 

17.1X10 4 

12.79 

3.33 

13.5 

4.04 


For all the data (see table 1) 


calc = 3. 6^6 Standard deviation, 0.32 


^Ergun^Ergun = 45^5 Standard deviation, 3.5 

For the minimum and maximum twisted-tape thickness and width, from table 2, and the combined data set of 
table 1, 



<w> 

(^exp) 

^calc/^Ergun 

Standard 

deviation 

0.318 cm (0.1252 in.) 

0.6276 

3.2 

0.28 


0.1296 cm (0.051 in.) 0.3287 cm (0.1294 in.) 0.5934 


NAS A/TM— 2002-2089 14 


29 





































































where (e exp ) is the area-averaged porosity and 


39 <JWm*Er S un<51 


However, the data sets for the Borda and orifice with screen and the orifice without screen have distinct flow 
characteristics, and more representative values of L|; nHm A r ,,. n)un can determined. For the 48-twisted-tape data of 
table 1, the coefficients of the least-squares trend line' 


lOglol^Ergunl -^1 lo g lol^ErguJ + !i \ 


through individual data sets are tabulated below, where R 2 is the regression coefficent. 


Flow configuration 

4, 

s, 

R 2 

Constant ~ ^Ergim -^Ergun 


Average porosity 

= 0.61 


Borda with screen 8 

-0.9873 

1.6673 

0.9965 

46.5 

Orifice with screen 

-1.0006 

1.6436 

.9984 

44.0 

Orifice without screen 

-0.9982 

1.6132 

.9988 

41.0 

Combined data sets 

-0.9867 

1.6469 

.992 

44.4 


Minimum porosity 

= 0.5934 


Borda with screen 8 

-0.9879 

1.6202 

0.9964 

41.7 

Orifice with screen 

-1.0006 

1.5958 

.9984 

39.4 

Orifice without screen 

-0.9982 

1.5655 

.9988 

36.8 

Combined data sets 

-0.9869 

1.5994 

.992 

39.8 


Maximum porosity 

= 0.6276 

Borda with screen 8 

-0.9866 

1.7190 

0.9966 

52.4 

Orifice with screen 

-1.0006 

1.6961 

.9984 

49.7 

Orifice without screen 

-0.9982 

1.6658 

.9988 

46.3 

Combined data sets 

-0.9864 

1.699 

.9919 

50.0 


8 One questionable data point set at average of previous and following points in table 1 . 


Estimates from figure 6 of Smithberg and Landis (1964) for air and water data are 




3 Re 


Ergun 


Exgun 


,= 3/ 
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Reynolds number. Re 

5000 

10 000 

20 000 

40 000 

60 000 

-^Ergun 

7500 

15 000 

30 000 

60 000 

90 000 

/at HJD = 3.62 

29x1 0 -3 

21xl0~ 3 

15.5xlO~ 3 

13xl0~ 3 

12.5xl0~ 3 

^Ergun 

87xl0~ 3 

63x10 3 

46.5xl0~ 3 

39x1 0 -3 

37.5xl0~ 3 

/at HJD = 4.34 

22.5x1 0 -3 

16.5x10 3 

13xl0~ 3 

1 lxl0~ 3 

9.7xl0~ 3 

^Ergun 

67.5xl0~ 3 

49.5xl0~ 3 

39x1 0 -3 

33xl0~ 3 

29-lxlO -3 

/at HJD = 10.3 

14xl0~ 3 

llxlO 3 

8.5X10 -3 

7.4X10’ 3 

6.8xl0~ 3 

^Ergun 

42x1 0 -3 

33x1 0 -3 

25.5xl0~ 3 

22.2xlO~ 3 

20.4X10" 3 

/at HJD = 22 

12xl0~ 3 

9.9x1 0 -3 

7.6x1 O’ 3 

6.6X10 -3 

6xl0" 3 

^Ergun 

36x1 O' 3 

29.7xl0~ 3 

22.8xl0~ 3 

19.8X10" 3 

18xl0 -3 

/ at HJD = oo 

8xlO -3 

7.1xl0' 3 

6xl0 -3 

5.35X10" 3 

5xlO~ 3 

^Ergun 

24x1 O' 3 

21.3xl0~ 3 

18xl0 -3 

16.1X10" 3 

15xlO -3 


The Smithberg and Landis (1964) data (N = 1 ) are Reynolds number dependent. It would be interesting to determine 
if data for N= 2, 3,..., 48 form parametric families that become more independent of Reynolds number as N in- 
creases (see eq. (A55». 

Estimates from Koch (1958) as provided from figure 7 of Smithberg and Landis (1964) are as follows: 



Reynolds number. Re 


2000 

3000 

6000 

10 000 

20 000 

30 000 

50 000 


^Ergun 


3000 

4500 

9000 

15 000 

30 000 

45 000 

75 000 

/at HID = 5 

32.5X10 -3 


18xl0~ 3 


11.6xl0“ 3 

lOxlO -3 

9xl0“ 3 

^Ergun 

97.5X10 -3 


54x1 0~ 3 


34.8X10 -3 


27xl0“ 3 

/at HID = 8.5 

21xl0 -3 


13.2X10 -3 

10.9X10 -3 

8.6xl0~ 3 

7.5xl0“ 3 

6.4xl0“ 3 

* Ergun 

63xl0“ 3 

52.5xl0~ 3 

39.6x1 0“ 3 

32.7x1 0 -3 

25.8X10 -3 

22.5x1 0“ 3 

19.2xl0“ 3 

fat HID =22 

14.2xl0~ 3 

12.3xl0~ 3 

9.7x1 0^ 3 

8.5X10' 3 

7.1X10 -3 

6.4xl0“ 3 

5.7xl0~ 3 

^Ergun 

42.6xl0~ 3 

36.9xl0~ 3 

29. lxlO 3 

25.5xl0~ 3 

21.3xl0~ 3 

19.2X10 -3 

17.1xl0“ 3 

f at HID — °° 

10xl0 -3 

9.5xl0~ 3 

8. Ixl0~ 3 

7.4xl0~ 3 

6.5xl0“ 3 

5.9xl0~ 3 

5.4X10' 3 

^Ergun 

30x1 0 -3 

28.5xl0~ 3 

24.3x1 0 -3 

22.2xl0~ 3 

19.5xl0~ 3 

17.7xl0~ 3 

16.2xl0“ 3 
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SUMMARY TABLE 



NASA/TM— 2002-2089 14 


32 














































NAS A/TM— 2002-2089 14 


33 


It is important to note that /Re = (function) = (conversion) ^Ergun^Ergun and Cj-= 4 f 0 . 
Turbulent calculations were approximately 30% lower than the data; equation not shown. 
Suggested general form is eq. (A85a). 

^eanalysis. 








APPENDIX B 


SYMBOLS 


A 


a 

a v 

Cf 

Co 

D 

D, 


D 


H 


D , 

D o 

D„ 


sphere 


tube 

(1 

(d) 


F(£) 


f 

fa 

fi 

fo 

f, 

G 


G 


o 


8 

80 

81 


8 2 

H 

H, 


K 

k 

k n 


flow area 


constant 

tube cross-sectional area without twisted tape 

surface area or area of sphere 

ratio of wetted surface area to bed volume 

specific surface area, x/o 

constant 

flow coefficient of friction 

constant relating velocity profiles to radial position 

tape diameter or width (for thin tubes, same as flow tube diameter) 

equivalent diameter 

hydraulic diameter; characteristic length of packed bed 
inside diameter 

virtual diameter of twisted tape 
equivalent particle diameter 
sphere diameter 
flow tube diameter 
fiber diameter 


average fiber diameter 

subscript denoting equivalent isothermal, eq. (A43) 
function relating friction to flow direction, eqs. (A97) 


( D Y dp / dx 

Fanning friction factor,^— 

no-swirl friction factor, eq. (A43b) 
friction factor, eq. (A75) 

Fanning friction factor for tube without twisted tape 


Darcy swirl (twisted tape) friction factor, f s = 4/ (Darcy/ = 4 Fanning/) 

mass flow, p u 

mass flux, p U a = WIA 0 

function defined in eq. (A24) 

function defined in eq. (A23) 

function defined in eqs. (A78) 

function defined in eq. (A82) 

tape twist through 1 80° or one-half full wave 

tape twist through 360° or one full wave 

torsion constant 

entrance pressure loss parameter, eq. (A47) 
packing constant, kJ'Q, eqs. (A86) and (A94) 
porosity shape parameter, eq. (A94) 
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k s surface roughness parameter 

L equivalent flow or twisted-tape length; nonideal flow path length 

L 0 tape length without twist; straight-line distance between pressure taps; bed or tube length 

M number much greater than unity 

m exponent, eq. (A42) 

m, number much less than unity 

N number of twisted tapes in tube 

n exponent, eqs. (A21) and (A43) 

n. number of 2 n twists 

P pressure, experimental 

p pressure, calculated 

R outer radius 

R h ratio of bed cross section available for flow to wetted perimeter 

Re Reynolds number, puD/ix 

Re equivalent Reynolds number 

Re ; Reynolds number, eq. (A75) 

r radial position 

r Q outer “wall” radius 

S wetted perimeter; surface area of fibrous specimen 

S 0 surface area, eq. (A101) 

? tape thickness 

(?) average twisted-tape thickness 

U bulk average axial velocity 

U e equivalent fluid velocity 

U 0 empty or unpacked-bed velocity 

u velocity, general or vector; superficial velocity 

M bar average velocity 

V bed volume 

V volumetric flow rate 

V s volume of sphere 

y sohd volume of solid in porous bed 

V total tota l bed volume (solid plus void) 

v flow velocity in fibrous specimen 

v 0 circumferential velocity; tangential fluid velocity 

v 0 velocity at r 0 

v t tangential velocity 

W mass flow rate 

w twisted-tape width 

(w) average twisted-tape width 

■^Ergun Ergun Reynolds number parameter, G 0 D p H 1 - e)j,t 

X n pressure drop data parameter 

x axial position 

7 ca]c calculated Ergun parameter 
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y Frgim Ergun friction factor parameter, (pA/ 5 / G„ |e 3 / (1 - e)| (p p / L 0 j 
Y n flow data parameter 

y number of tube diameters per 180° of twist (also used as coordinate, fig. 2(e)) 

Z analytical axial locus 

Z g equivalent axial distance without twist 

z axial coordinate 

a, <]), % polar coordinate parameters, fig. 2(e) 

8 surface roughness 

j 0 bulk density parameter, eq. (A86) 

e bed porosity parameter 

e experimental bed porosity 

E m odei modeled bed porosity 

k torsion parameter, eq. (A71) 

p average fluid density 

p viscosity 

P 

v ~ kinematic velocity 

(p analytical twist parameter 

(|) angle between interface normal and microscopic flow 

x/o surf ace- to- volume parameter, eq. (A86) 

x thickness of fibrous specimen 

a deviation 

C, porosity parameter, eq. (A94) 
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TABLE 1.— FLOW AND PRESSURE DROP DATA AND ANALYSIS FOR 48 TWISTED TAPES IN CYLINDRICAL TUBE 


V, 

cm 3 /min 

A P, 
psi 

G 0 . , 

g/cm 2 -s 

dP, 

g/cm-s 2 

^Ergun 

0.659 G 0 

^Ergun 

73.34x10^ 
dPIG„ 2 

y — 

1 calc 

150/X^ UII + 

-^calcAErgun 


Y Y 

^-Ergun 1 Ergun 

0 2 

a 8218 

3.55 

24.33 

244762 

15.28 

1.81 

11.57 

6.38 

0.0039 

27.69 

3.42 

7657 

3.29 

22.66 

226836 

14.23 

1.94 

12.29 

6.35 

0.0009 

27.54 

2.89 

5900 

2.54 

17.46 

175125 

10.97 

2.52 

15.43 

6.13 

0.0356 

27.60 

3.08 

4730 

2.07 

14.00 

142720 

8.79 

3.19 

18.81 

5.90 

0.1799 

28.05 

4.89 

3881 

1.67 

11.49 

115141 

7.21 

3.82 

22.54 

5.90 

0.1797 

27.58 

3.03 

2804 

1.11 

8.30 

76531 

5.21 

4.87 

30.53 

6.27 

0.0024 

25.37 

0.22 

2790 

1.14 

8.26 

78600 

5.19 

5.05 

30.67 

6.07 

0.0606 

26.19 

0.12 

1546 

0.59 

4.58 

40679 

2.87 

8.51 

53.95 

6.34 

0.0003 

24.46 

1.90 

1272 

0.47 

3.77 

32405 

2.36 

10.02 

65.19 

6.51 

0.0353 

23.68 

4.65 

b 957 

0.31 

2.83 

21374 

1.78 

11.67 

86.07 

7.37 

1.1111 

20.76 

25.77 

975 

0.37 

2.89 

25510 

1.81 

13.42 

84.51 

6.30 

0.0005 

24.32 

2.30 

1518 

0.63 

4.49 

43437 

2.82 

9.43 

54.91 

5.82 

0.2458 

26.60 

0.58 

2014 

0.87 

5.96 

59984 

3.74 

7.40 

41.82 

5.65 

0.4437 

27.69 

3.42 

2713 

1.13 

8.03 

77910 

5.04 

5.29 

31.49 

5.95 

0.1377 

26.70 

0.74 

3538 

1.59 

10.47 

109626 

6.58 

4.38 

24.56 

5.61 

0.5088 

28.81 

8.80 

4307 

1.94 

12.75 

133757 

8.01 

3.61 

20.49 

5.68 

0.4088 

28.87 

9.19 

5260 

2.33 

15.57 

160647 

9.78 

2.90 

17.09 

5.89 

0.1888 

28.39 

6.52 

6129 

2.67 

18.14 

184088 

11.39 

2.45 

14.92 

6.09 

0.0548 

27.92 

4.34 

7151 

3.09 

21.17 

213046 

13.29 

2.08 

13.03 

6.26 

0.0041 

27.70 

3.45 

8207 

3.49 

24.29 

240625 

15.26 

1.79 

11.58 

6.48 

0.0264 

27.26 

2.01 

6354 

2.73 

18.81 

188225 

11.81 

2.33 

14.45 

6.20 

0.0151 

27.54 

2.89 

4435 

1.92 

13.13 

132378 

8.24 

3.37 

19.94 

5.93 

0.1557 

27.75 

3.65 

2471 

1.02 

7.31 

70326 

4.59 

5.76 

34.41 

5.97 

0.1205 

26.46 

0.38 

c 1592 

0.686 

4.71 

47298 

2.96 

9.33 

52.44 

5.62 

0.4925 

27.62 

3.17 

2721 

1.174 

8.05 

80944 

5.06 

5.47 

31.41 

5.74 

0.3320 

27.66 

3.30 

4349 

1.866 

12.87 

128655 

8.08 

3.40 

20.30 

5.97 

0.1236 

27.50 

2.76 

4311 

1.84 

12.76 

126862 

8.01 

3.41 

20.47 

6.00 

0.1054 

27.36 

2.31 

6112 

2.58 

18.09 

177883 

11.36 

2.38 

14.95 

6.28 

0.0017 

27.06 

1.48 

8045 

3.37 

23.81 

232351 

14.95 

1.80 

11.78 

6.56 

0.0581 

26.85 

1.02 

7193 

2.96 

21.29 

204083 

13.37 

1.97 

12.97 

6.57 

0.0644 

26.38 

0.29 

6057 

2.49 

17.93 

171678 

11.26 

2.34 

15.07 

6.44 

0.0144 

26.35 

0.26 

4993 

2.04 

14.78 

140652 

9.28 

2.82 

17.91 

6.35 

0.0008 

26.19 

0.12 

3459 

1.41 

10.24 

97215 

6.43 

4.06 

25.08 

6.17 

0.0221 

26.13 

0.08 

2298 

1.04 

6.80 

71705 

4.27 

6.79 

36.86 

5.43 

0.7948 

29.01 

10.04 

1526 

0.64 

4.52 

44126 

2.84 

9.48 

54.63 

5.76 

0.3087 

26.88 

1.09 

958 

0.39 

2.84 

26889 

1.78 

14.65 

85.98 

5.87 

0.2045 

26.09 

0.06 

446 

0.2 

1.32 

13789 

0.83 

34.67 

182.68 

5.27 

1.1047 

28.74 

8.43 

d 5859 

2.52 

17.34 

173746 

10.89 

2.53 

15.52 

6.13 

0.0353 

27.57 

2.99 

8156 

3.48 

24.14 

239936 

15.16 

1.80 

11.64 

6.45 

0.0181 

27.35 

2.28 

6777 

2.86 

20.06 

197188 

12.60 

2.15 

13.66 

6.36 

0.0016 

27.05 

1.47 


a Borda + screen (12/10/98). 

Questionable. 

c Borda + screen (12/12/98). 

d Borda + screen (12/14/98). 

e Orifice + screen (12/14/98; reverse of Borda + screen) . 
f Orifice without screen (12/14/98; reverse of Borda without screen). 
8 Orifice without screen (12/19/98; reverse of Borda without screen). 
h N - 1 points. 

W points. 

J Throw out one point. 
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TABLE L— CONCLUDED. 



A P, 
psi 

G 0 , 

g/cm 2 -s 

dP, 

g/cm-s 2 

-^Ergun — 

0.659 G 0 

73.34x10^ 

dPIG* 

y^ajc = 

150/A Ergu „+ 

1.75 

^calc/^Ergun 


v v 

^-Ergun 1 Ergun 

a 2 

4546 

1.9 

13.46 

130999 

8.45 

3.17 

19.50 

6.15 

0.0285 

26.79 

0.90 

2367 

1 

7.01 

68947 

4.40 

6.15 

35.84 

5.82 

0.2466 

27.08 

1.54 

1768 

0.75 

5.23 

51710 

3.29 

8.27 

47.39 

5.73 

0.3505 

27.19 

1.83 

1245 

0.55 

3.69 

37921 

2.31 

12.24 

66.56 

5.44 

0.7740 

28.32 

6.13 

833 

0.36 

2.47 

24821 

1.55 

17.89 

98.62 

5.51 

0.6519 

27.70 

3.47 

718 

0.32 

2.13 

22063 

1.33 

21.40 

114.14 

5.33 

0.9753 

28.57 

7.44 

e 982 

0.35 

2.91 

24131 

1.83 

12.52 

83.92 

6.71 

0.1487 

22.85 

8.97 

2239 

0.88 

6.63 

60673 

4.16 

6.05 

37.79 

6.24 

0.0059 

25.19 

0.42 

4754 

1.88 

14.07 

129620 

8.84 

2.87 

18.72 

6.53 

0.0431 

25.35 

0.24 

8239 

3.27 

24.39 

225457 

15.32 

1.66 

11.54 

6.95 

0.3964 

25.44 

0.16 

7030 

2.76 

20.81 

190294 

13.07 

1.93 

13.23 

6.87 

0.3017 

25.17 

0.46 

4265 

1.65 

12.62 

113763 

7.93 

3.13 

20.67 

6.61 

0.0832 

24.80 

1.09 

2000 

0.77 

5.92 

53089 

3.72 

6.64 

42.10 

6.34 

0.0005 

24.68 

1.35 

1072 

0.42 

3.17 

28958 

1.99 

12.60 

77.02 

6.11 

0.0434 

25.11 

0.53 

722 

0.28 

2.14 

19305 

1.34 

18.52 

113.51 

6.13 

0.0366 

24.86 

0.96 

452 

f 1297 

0.19 

1.34 

13100 

0.84 

32.07 

180.28 

5.62 

0.4877 

26.94 

1.22 

0.47 

3.84 

32405 

2.41 

9.63 

63.97 

6.64 

0.1020 

23.23 

6.82 

3303 

1.19 

9.78 

82047 

6.14 

3.76 

26.18 

6.96 

0.4104 

23.09 

7.54 

5038 

1.83 

14.91 

126173 

9.37 

2.49 

17.77 

7.15 

0.6827 

23.28 

6.54 

6767 

2.47 

20.03 

170299 

12.58 

1.86 

13.67 

7.35 

1.0651 

23.40 

5.97 

8412 

3.08 

24.90 

212357 

15.64 

1.50 

11.34 

7.56 

1.5308 

23.47 

5.62 

7619 

2.76 

22.55 

190294 

14.16 

1.64 

12.34 

7.53 

1.4576 

23.22 

6.87 

6019 

2.15 

17.82 

148236 

11.19 

2.05 

15.16 

7.41 

1.1803 

22.90 

8.67 

4100 

1.45 

12.14 

99973 

7.62 

2.97 

21.43 

7.21 

0.7837 

22.67 

10.06 

2003 

0.704 

5.93 

48539 

3.72 

6.05 

42.04 

6.95 

0.3935 

22.53 

10.96 

1002 

0.34 

2.97 

23442 

1.86 

1 1.68 

82.28 

7.05 

0.5276 

21.75 

16.73 

868 

0.31 

2.57 

21374 

1.61 

14.19 

94.72 

6.68 

0.1265 

22.89 

8.69 

556 

0.208 

1.65 

14341 

1.03 

23.20 

1 46.SS 

6.33 

0.0001 

23.98 

3.46 

8 10 13 

0.389 

3.00 

26820 

1.88 

13.07 

81.41 

6.23 

0.0085 

24.61 

1.50 

3095 

1.17 

9.16 

80668 

5.75 

4.21 

27.82 

6.61 

0.0817 

24.23 

2.59 

6556 

2.485 

19.41 

171333 

12.19 

1.99 

14.06 

7.05 

0.5352 

24.30 

2.38 

8239 

3.123 

24.39 

215321 

15.32 

1.59 

11.54 

7.28 

0.9154 

24.30 

2.38 

5000 

1.851 

14.80 

127621 

9.29 

2.55 

17.89 

7.01 

0.4715 

23.73 

4.46 

1700 

0.621 

5.03 

42816 

3.16 

7.41 

49.22 

6.64 

0.1039 

23.42 

5.88 

743 

0.281 

2.20 

19374 

1.38 

17.55 

110.36 

6.29 

0.0011 

24.24 

2.55 






#N/A 


474.00 

22.56 

1938.03 

295.77 








6.32 

0.3048 

25.84 

h 4.00 









0.3007 


*3.94 

J 3. 648663 


a Borda + screen (12/10/98). 

b Questionable. 

c Borda + screen (12/12/98). 

d Borda + screen (12/14/98). 

e Orifice + screen (12/14/98; reverse of Borda + screen) . 
f Orifice without screen (12/14/98; reverse of Borda without screen). 
8 Orifice without screen (12/19/98; reverse of Borda without screen). 
h N - 1 points. 
l N points. 

J Throw out one point. 
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TABLE 2.— TWISTED-TAPE WIDTH <w> AND THICKNESS ( t ) MEASUREMENTS 
FOR 48 TWISTED TAPES a 
[Sorted in ascending order (estimated).] 


Tape 

Left-end measurements 

Right- end measurements 


Thickness, 

Width, 

A= wt. 

A= wt, 

Thickness, 

Width, 



M, 

. 2 

in. 

. 2 

in. 

<0. 

M. 



in. 

in. 



in. 

in. 

1 

0.038 

0.129 

0.004902 

0.004514 

0.037 

0.122 

2 

.038 

.124 

.004712 

.004674 

.038 

.123 

3 

.038 

.123 

.004674 

.004636 

.038 

.122 

4 

.039 

.122 

.004758 

.004674 

.038 

.123 

5 

.039 

.123 

.004797 

.005265 

.039 

.135 

6 

.039 

.125 

.004875 

.004797 

.039 

.123 

7 

.04 

.129 

.00516 

.00488 

.04 

.122 

8 

.04 

.124 

.00496 

.00492 

.04 

.123 

9 

.041 

.126 

.005166 

.005002 

.041 

.122 

10 

.042 

.118 

.004956 

.005453 

.041 

.133 

11 



.124 

.005208 

.00504 

.042 

.12 

12 



.118 

.004956 

.005418 



.129 

13 



.129 

.005418 

.005292 



.126 

14 



.128 

.005376 

.005208 



.124 

15 



.124 

.005208 

.004914 

1 

r 

.117 

16 

1 

r 

.128 

.005376 

.005504 

.043 

.128 

17 

.043 

.121 

.005203 

.005192 

.044 

.118 

18 

.043 

.121 

.005203 

.00528 

.044 

.12 

19 

.043 

.119 

.005117 

.005324 

.044 

.121 

20 

.044 

.119 

.005236 

.005715 

.045 

.127 

21 

.047 

.123 

.005781 

.00611 

.047 

.13 

22 

.048 

.145 

.00696 

.0066 

.05 

.132 

23 

.049 

.13 

.00637 

.0067 

.05 

.134 

24 

.05 

.124 

.0062 

.0064 

.05 

.128 

25 

.051 

.125 

.006375 

.006528 

.051 

.128 

26 

.052 

.124 

.006448 

.0065 

.052 

.125 

27 



.128 

.006656 

.006916 



.133 

28 



.123 

.006396 

.006552 



.126 

29 



.126 

.006552 

.0065 



.125 

30 



.125 

.0065 

.007182 

3 

f 

.133 

31 



.125 

.0065 

.00767 

.059 

.13 

32 



.124 

.006448 

.00944 

.059 

.16 

33 


r 

.125 

.0065 

.007847 

.059 

.133 

34 

.053 

.126 

.006678 

.008296 

.061 

.136 

35 

.053 

.125 

.006625 

.008235 



.135 

36 

.054 

.127 

.006858 

.007503 



.123 

37 



.123 

.006642 

.007747 

1 

f 

.127 

38 



.126 

.006804 

.007936 

.062 

.128 

39 



.123 

.006642 

.008618 

.062 

.139 

40 

1 

r 

.125 

.00675 

.008866 

.062 

.143 


a Average porosity, 0.6105. Corresponding average t (in.) = 0.0502 and w (in.) = 0.1273. 
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TABLE 2.— CONCLUDED. 


Tape 

Left-end measurements 

Right-end measurements 


Thickness, 

Width, 

A= wt. 

A= wt. 

Thickness, 

Width, 


(f). 

M, 

. 2 

in. 

. 2 

in. 

« >. 

M, 


in. 

in. 



in. 

in. 

41 

0.055 

0.132 

0.00726 

0.008253 

0.063 

0.131 

42 

.055 

.126 

.00693 

.008064 

.063 

.128 

43 

.057 

.125 

.007125 

.009009 

.063 

.143 

44 

.061 

.127 

.007747 

.00896 

.064 

.14 

45 

.061 

.132 

.008052 

.00896 

.064 

.14 

46 

.062 

.118 

.007316 

.008704 

.064 

.136 

47 

.062 

.126 

.007812 

.00871 

.065 

.134 

48 

.065 

.128 

.00832 

.008844 

.066 

.134 

Totals 

2.334 

6.01 

0.292508 

0.319352 

2.45 

6.212 

Average 

0.04862 

0.12521 


0.05104 

0.12942 

Porosity 


0.6276 

0.5934 



a Average porosity, 0.6105. Corresponding average t (in.) = 0.0502 and w (in.) = 0.1273. 



Figure 1 . — Single twisted tapes. (From Smithberg and 
Landis, 1 964.) 
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Figure 2. — Configurations of packed bed of twisted tapes, (a) Typical twisted-tape assembly, (b) Cross 
section (end view), (c) Ideal packed bed in soccer ball configuration, (d) Test section, (e) Test section 
screen, (f) Twisted tapes, (g) Prismatic bar model. 
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Figure 3. — Behavior for single twisted tape and 48 twisted tapes in packed bed relative to 
Ergun model for laminar and turbulent-flow data. 


NAS A/TM— 2002-2089 14 


47 




' Ergun 







NAS A/TM— 2002-2089 14 


49 


Figure 6. — Axial flow field along centerline of packed bed of twisted tapes, (a) Inlet region, (b) Developed region. 


Ergun friction factor, VE r g Un 



Figure 7. — Transverse flow field at one diameter from packed-bed inlet. 



1CT 1 10° 10 1 1 


Ergun Reynolds number, X Ergun 

Figure 8. — Laminar flow behavior for single twisted tape and packed bed of 48 twisted 
tapes. 
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Ergun friction factor, / Ergun Ergun friction factor, y E| g Un 



Ergun Reynolds number, X Ergun 

Figure 9. — Velocity-corrected turbulent data and correlation functions. 
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10” 1 


10-2 
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Ergun Reynolds number, X Ergun 

Figure 1 0. — Behavior of single twisted tape and packed bed of 48 twisted tapes relative to 
Ergun model for laminar flow data. 
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Fiber concentration fraction, 


° uytr* + a8 = k ° (eq - (A95)) 

O Glass wool 

® Glass wool and copper wire 
Cl Glass (fibers perpendicular to flow) 
£1 Glass (fibers parallel to flow) 

□ Cotton wool 



\Pd z /fxLU 0 

Figure 12. — Modified fit of Davies data (1952). (From fig. 3 of Hersh and Walker, 

1980) APd 2 ~| „ 1 6 k 0 (1 - A 2 APd 2 

fxLU° J para « e | ~ « 


normal 


NAS A/TM— 2002-2089 14 


53 



^SlO^Ergun -n lo 9lO y Ergun 



10 ° 10 1 10 2 10 3 10 4 10 5 

lo 9io X Ergun 

igure 13. — Dependence of Ergun model including packed-sphere, turbulent-flow data of 
Wentz and Thodos (1 963). 
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□ 0 


HI 


- 2.0 


0.03(l50/X Ergun + 4.5/X^ n )l Thodos (1963) 


-2.5 


□ X Ergun /(l/ e /U), Smith berg and Landis (1964); Koch (1 958) 
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Figure 1 4. — Relative magnitudes and slopes of flow parameters for packed beds of spheres 
(Wentz and Thodos, 1 963) compared with data for single twisted tape in tube (Smithberg 
and Landis, 1 964, and Koch, 1 958). 
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